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Abstract 

We calculate the probability that the rapidity gaps in diffractive processes survive both 
eikonal and enhanced rescattering. We present arguments that enhanced rescattering, 
which violates soft-hard factorization, is not very strong. Accounting for NLO effects, 
there is no reason to expect that the black disc regime is reached at the LHC. We discuss 
the predictions for the survival of the rapidity gaps for exclusive Higgs production at the 
LHC. 



1 Overview of rapidity gap survival 

We discuss the effect of screening corrections (that is, the rapidity gap survival factors) in 
diffractive processes. In particular we study the exclusive central production of Higgs bosons 
or another small-size heavy system, which we collectively denote by A. The process we have in 
mind is pp p + A + p, where the + signs indicate rapidity gaps. Usually the cross sections of 
such processes are calculated assuming soft-hard factorization. Symbolically it can be written 
as 

^ = \fg^MA®fg\^-S^ (1) 




Figure 1: The mechanism for the exclusive process pp p + A + p, with the eikonal and 
enhanced survival factors shown symbolically. The thick lines on the Pomeron ladders, either 
side of the subprocess [gg — > A), indicate the rapidity interval Ay where enhanced absorption 
is not permitted [2, 3, 4]. 

where the first factor is the bare amplitude for the exclusive process calculated in perturbative 
QCD as the convolution of the hard matrix element Ma and the gluon distributions fg. This 
is shown in the right half of Fig. 1. The factor S"^ is the probability that the secondaries, which 
are produced by soft rescattering, do not populate the gap. 

As mentioned above, the cross section in form (1) assumes soft-hard factorization. The 
survival factor, calculated from a model of soft interactions, does not depend on the structure 
of the original perturbative QCD amplitude; symbolically it is shown by Seik in the left part of 
the figure. 

Actually, the situation is more complicated. We need to consider so-called "enhanced" 
rescattering corrections which involve intermediate partons, as indicated by the typical contri- 
bution to Sgjjj^ shown in Fig. 1. As a rule, the partons which participate in the hard subprocess 
have large virtuality in transverse momentum and the probability of rescattering of such partons 
is negligible. Due to the strong /cj-ordering, the most important rescattering involves partons 
near the beginning of the evolution, that is near the input scale go- However, the value of go is 
not fixed. To be able to neglect all the enhanced absorptive corrections we have to choose large 
go. On the other hand, the development of the parton shower starts from a low scale gg ~ l/R, 
where R is the proton radius. This part of the evolution of the parton shower from the beam 
may be affected by the interaction with the target (and vice versa). In this way 'enhanced' 
rescattering violates soft-hard factorization, leading to a distortion of the input distribution at 
go. Our aim is to quantify the rescattering effect which originates in the interval [q^, go). To 
do this we need to introduce components, a, of the Pomeron of different size. Here, we take 
three components, a = Pi, P2, P3 corresponding to large-, intermediate- and small-size as in 
the preceding paper [1]. These were associated with transverse momenta of about 0.5, 1.5 and 
5 GeV respectively. 

The opacity of the 'target A;' (or beam i) proton with respect to intermediate parton c, 
generated by Pomeron component a is denoted by XQ^ (or XQ^). The parameter A relates the 
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Figure 2: (a) The ladder structure of the triple-Pomeron amphtude between diffractive eigen- 
states of the proton; the rapidity y spans an interval to y = Ins. (b) A multi- 

Pomeron diagram. 



triple-Pomeron coupling to the Pomeron coupling to the proton: gppp = Xg^- Its presence 
will allow multi- Pomeron- Pomeron absorptive effects to be included [1]. If, for the moment, we 
consider a Pomeron with only one component, then the opacity satisfies the evolution equation 
in rapidity y, at impact parameter b, 

dnjy,b) ^ ^.xn.iy'M)/2 ^-xnuy,m + a'^ j 6) , (2) 

where y' = Ins— Let us explain the meanings of the three factors on the right-hand-side of (2). 
If only the last factor, {...)Qk, is present then the evolution generates the ladder-type structure 
of the bare Pomeron exchange amplitude, where the Pomeron trajectory ap = 1 + A + a't. 
The inclusion of the preceding factor allows for rescatterings of an intermediate parton c with 
the "target" proton k; Fig. 2(a) shows the simplest (single) rescattering which generates the 
triple-Pomeron diagram. Finally, the first factor allows for rescatterings with the beam i. In 
this way the absorptive effects generated by all multi- Pomeron- Pomeron diagrams are included, 
like the one shown in Fig. 2(b) containing, in general, m — > n multi-Pomeron vertices. There 
is an analogous equation for the evolution of Qi{y',b), and the two equations may be solved 
iteratively. 

Now, consider a Pomeron with more components a, then these equations have a matrix 
form in aa' space [1]. At each step of the evolution the equation includes absorptive factors 
of the form e"'^^"*^^"''^/^. By solving the equations with and without these suppression factors, 
we can quantify the effect of enhanced absorption. Moreover, we start the evolution from the 
large component Pi, and since the evolution equations allow a transition from one component 
to another (corresponding to BFKL diffusion [5] in Inkt space), we can see how the enhanced 
absorption will affect the distribution of high kt in the small size component P3. 

If we consider pure elastic eikonal rescattering. Fig. 3, then phenomeno logically we can 
describe the bare high-energy elastic amplitude F by just one Regge pole, see Fig. 4(a). On the 
other hand, we know that a sizeable part of the cross section comes from diffractive dissociation. 
High-mass, M, diffractive dissociation is usually described by triple-Pomeron diagram, and we 
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Figure 3: Pure elastic eikonal scattering. 



(a) 



(b) 



(c) 




M 



F = 



+ 




A 



Figure 4: (a) The bare pole amplitude; (b) the triple-Pomeron diagram describing high mass, M, 
diffractive dissociation; (c) the contribution of the triple-Pomeron coupling to central exclusive 
production. 

have to include this contribution both in the irreducible amplitude F (that is, the pp opacity) 
as in Fig. 4(b), and in the bare amplitude for central exclusive production, as in Fig. 4(c). The 
latter corresponds to a contribution to the so-called enhanced screening factor S^nh of Fig 1. 
Physically, this factor accounts for the absorption of intermediate partons inside the original 
central exclusive amplitude. It is called enhanced since we have to multiply the probabilities of 
absorption on each individual intermediate parton, and thus the final effect is enhanced by the 
large multiplicity of intermediate partons. 

A complication is that Scnh cannot be considered simply as an overall multiplicative factor. 
The probability of interaction with a given intermediate parton depends on its position in 
configuration space; that is, on its impact parameter b and its momentum fcj. This means that 
S'enh simultaneously changes the distribution of the active partons which finally participate in 
the hard subprocess. 

We are now in a position to calculate the rapidity gap survival factors including not just 
eikonal rescattering, but also allowing for the possibility of enhanced screening. 
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2 Eikonal screening 



The gap survival factor caused by eikonal rescattering of the Good- Walker eigenstates [6] , for 
a fixed impact paramter b, is 



i,k 


2 

\ak\'^ Mik(h) exp(-Fifc(s,b)/2) 




X]|ai|2 |ajk|2 Mik{h) 

i,k 


2 



(3) 



The Fik are the s-channel two-particle irreducible amplitudes for the scattering of the various 
diffractive eigenstates^ and (pk), for given separations b = bi — b2 between the incoming 
protons 



(4) 



where Y = Ins. It is the total opacity of the ik interaction. It includes the sum over the 
different Pomeron components a and the integral over the impact parameter which accounts 
for the overlap of the partonic wave functions (that is, the opacities and flk) generated by 
the beam and target protons for fixed b = bi — b2 separation. The factor 1/Pq provides the 
appropriate normalisation as explained in the preceding paper [1] . 

The most topical case is Central Exclusive Diffractive (CED) production of a Higgs boson 
(or other heavy system of small size) [7, 8, 9, 10]. Then the matrix element Mik of the hard 
subprocess gg ^ A can be calculated perturbatively, since the large mass and small size of the 
centrally produced object introduces a large scale [7, 8]. However, the shape of Adik in impact 
parameter space b and the relative couplings of the hard matrix element to different Good- 
Walker eigenstates are not known. Their calculation is beyond the ability of perturbative 
QCD. 

One possibility is to say that the b dependence of M. should be, more or less, the same as 
for diffractive J /ip electroproduction (7+p J/ip+p), and the coupling to the component 
of the proton should be proportional to the same factor 7^, as in a soft interaction. This leads 
to 

Mik oc 7i7fc exp(-6V45) (5) 

with 5 ~ 4 GcV^^ [11]. Under these assumptions, the 'soft' model of [1] predicts S^^]^=0.017 
for the LHC energy, which is close to the value 0.023 obtained recently by the Tel-Aviv group 
[12]. 



^Recall that, in the global model of soft interactions of the preceding paper [1], the Good- Walker formalism 
[6] was used to decompose the proton wave function in terms of three diffractive eigenstates \p) = J2i^i\'f'i)- 
That is, we used a three-channel eikonal to describe elastic rescattering and low-mass proton dissociation. 
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Figure 5: (a): A symbolic representation of (6); (b) and (c): diagrams relevant for the compu- 
tation of Sgj^i^ and S^^^S^^^ of (8) and (9), where is the opacity when the screening corrections 
embracing the hard matrix element are neglected. 

Expressing the survival factors in this manner is too simplistic. However, these numbers are 
frequently used as the reference point. This can be misleading, since S"^ calculated in this way 
strongly depends on an additional assumption about the slope B, a quantity which is mainly 
of 'soft' origin. 

Another, more plausible and self-consistent, possibility is to calculate A^jfc(b) under the 
assumption that the active gluons that participate in the hard subprocess come from the small- 
size component P3 of the Pomeron. That is, 

A^^f (b) = V / n[%yiMM)n^'{y2MM)d\d%6^'\h^-h,-h) , (6) 

where the rapidities yj = — In^j with j = 1,2, correspond to the momentum fractions carried by 
the active gluons which participate in the hard subprocess with matrix element M{gg A). 
The superscript "enh" is to indicate that the matrix element calculated in this way already 
accounts for the "enhanced" rescattering on intermediate partons. A symbolic representation 
of the formula is given in Fig. 5(a). The formula is similar to formula (4) for the irreducible 
amplitude Fjfc(b). The difference is that now the two active gluons are separated by a large 
rapidity interval InM^, so that yi + y2 = Y — InM^. 

3 Enhanced screening 

However, as mentioned above, besides pure eikonal screening, there may be enhanced rescat- 
tering arising from absorption of the intermediate partons in the hard matrix element Aitk- 
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Working at LO (of coUinear approximation) we have to neglect such an effect. Due to strong 
/CfOrdering the transverse momenta of all the intermediate partons are very large (i.e. the 
transverse size of the Pomeron is very small) and therefore the absorptive effects are negligible. 
Nevertheless, this may be not true at a very low x, where the parton densities become close 
to saturation and the small value of the absorptive cross section is compensated by the large 
value of the parton density. 

Indeed, the contribution of the first enhanced diagram, which describes the absorption of an 
intermediate parton, was estimated in the framework of the pcrturbative QCD in Ref.[13]. It 
turns out that it could be quite large. On the other hand, such an effect is not seen experimen- 
tally. The absorptive corrections due to enhanced screening must increase with energy. This 
is not observed in the present data (see [3] for a more detail discussion). One reason is that 
the gap survival factor S^^^ already absorbs almost all the contribution from the centre of the 
disk. The parton essentially only survives eikonal rescattering on the periphery, that is at large 
impact parameters h. On the other hand, on the periphery, the parton density is rather small 
and the probability of enhanced absorption is not large. This fact can be seen in Ref. [14]. 
There the momentum, Q^, below which we may approach saturation, was extracted from the 
HERA data in the framework of the dipole model. Already at 6 = 0.6 fm the value of Q\ < 0.3 
GeV^ for the relevant x ~ 10~^. 

Recall that, in our model, the absorption caused by rescattering of intermediate partons 
is included as exponential factors, Q-K^i+^k)/"^^ jj-^ equations (2) for the opacities. So we have 
the possibility to quantify the role of this ('enhanced') exponential suppression by solving the 
evolution equations with and without the exponential factors. 

First, we could calculate the value of A4ik excluding any effect of 'enhanced' absorption. 
To do this we would delete the factors e~-^(^'=+^')/^ from the equation (2), and then solve the 
modified equations for opacities ilj and respectively. However, this is not what we need, 
and before we start the calculation it is necessary to note two pragmatic observations. The 
problem is that we no longer have exact factorisation between the hard and soft parts of the 
process. Thus before computing the effect of soft absorption wc must fix what is included in 
the bare CED amplitude calculated in terms of pcrturbative QCD. 

4 Matching 5^ with the hard matrix element 

The first observation is that the bare amplitude is calculated as a convolution of two gener- 
alised (skewed) gluon distributions with the hard subprocess matrix element [7, 8] . These gluon 
distributions are determined from integrated gluon distributions of a global parton analysis of 
mainly deep inelastic scattering data. However, these phenomenological integrated parton dis- 
tributions already include the interactions of the intermediate partons with the parent proton. 
As a consequence, in order to determine S'^^jj we should compare the matrix element calculated 
including all enhanced diagrams, like the one sketched in Fig. 5(a), with the matrix element 
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obtained by excluding, in the evolution of fij, the exponential factor e~'^^'=/^ which 'embraces' 
the hard STibproccss, and keep only e"'^^*/^, as indicated in Fig. 5(b); and vice versa for the 
evolution of Qk- In this way we mimic the gluons used in the perturbativc QCD calculation of 
the cross section. Let us denote the opacities, obtained this way, as Q and the corresponding 
matrix element of (6) as Ai,;f..(h). These opacities mimic the phenomenological gluon distribu- 
tions used in the perturbativc QCD calculation. So, for fixed indices i, k, the enhanced survival 
factor is given by 

^enh 



c;r2 
'-'enh 



M 



(7) 



Now we are ready to calculate the survival factors. Wc would like to compare the survival 
factors computed as indicated in Fig. 5(b) and (c). First, corresponding to Fig. 5(c), we 
compute 5'^ji^(b) from (3) where A^^^'^(b) includes the full enhanced screening, but not the 
eikonal rescattering. Symbolically we have 



M 



enh 



* e 



-Fik/2 



M 



enh 



(8) 



where TVf'^'^'^' is given by (6). Then wc compute the total absorptive effect (Sgj^^C^) 'S'gi]^(b) by 
replacing M'^"^ in the denominator by A4. The calculation of corresponds to Fig. 5(b). 
That is. 



M 



enh 



* e 



-F^k/2 



M 



(9) 



The second observation is that the phenomenologically determined generalised gluon distri- 
butions are usually taken at pt — 0, and then the observed "total" cross section is calculated 
by integrating over pt of the recoil protons assuming an exponential behaviour e~^^* ; that is 



/ 



da 



dpudp 



-dpltdplt = 



da 



'2t 



dplMt 



Plt=P2t=0 



where 



/ dp', e-^^? = 1/B = (p?). 



(10) 



(11) 



However, the total soft absorptive effect changes the pt distribution in comparison to that for 
the bare cross section determined from perturbativc QCD. Moreover, the correct pt dependence 
of matrix element (6) does not have an exponential form. Thus the additional factor introduced 
by the soft interactions is not just the gap survival S'^, but rather S'{p,)'', where the square 
arises since we have to integrate over the pt distributions of two outgoing protons. Indeed in 
all the previous calculations the soft pref actor had the form^ S'^/B'^. Now, in our model the 



■^At larger impact parameter b the absorption is weaker. Hence the value of increases with the slope B. 
It was shown that the ratio S^/B^ is approximately stable for reasonable variations of B [15]. 
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Figure 6: The survival factors for the Central Exclusive DifFractive (CED) production of a 
heavy 0+ system of mass M^=150 GeV and Ma= 10 GeV. 



b behaviour of M.{h) is driven by the opacities. Thus we present the final result in the form 
S'^{pf)'^, and so to compare it with previous predictions obtained for S = 4 GeV~^ we need to 
introduce the "renormalisation" factor {{pf)BY. The resulting (effective) value is denoted by 



5 Impact parameter dependence of S"^^-^ 

First, we present the results of the calculation of the survival factors at fixed pp impact param- 
eter b = bi — that is, before integrating over d'^b. The results are presented in Fig. 6 for 
the LHC energy y/s = 14 TeV for, first, the production of a Higgs boson of mass Ma = 120 
GeV and, second, for the central producton of an object of mass = 10 GeV (which may be 
a pair of high-ET photons or a Xb)- For the 10 GeV object we assume that the active gluons 
come from the intermediate Pomeron state, P2. As expected the eikonal suppression, S'gjj,, is 
almost the same for both objects (with a bit stronger suppression for P3 which lies in a bit 
smaller b domain). On the other hand, the enhanced suppression, S^^^i^ is stronger for the 
lighter object as the rapidity interval allowed for the interaction with the intermediate partons 
is considerably larger, and the absorption cross section is larger for a lower kt parton. 

After we multiply the bare QCD amplitude, shown by the dashed line in the plots of Fig. 6, 
by the total survival factors we obtain the impact parameter profiles of the final CED signal 
which peak around 0.8 fm. The signal vanishes as 6 — simply because we plot the integrand 
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6|Alp of / d'^b S'^\M(h)\'^ = Stt / b\M(h)\'^db; the normalisation of the signal in the figure was 
chosen to facilitate the plot. 

Assuming that all the partons are distributed homogeneously, we would expect S^^-^ — > 1 
at the large values of b of the periphery. However, the results show only a weak tendency to 
increase at large b. Thus, in the 'soft' model of [1] (with a small value of a'p), S^^jj comes mainly 
from "hot-spots" in which many individual intermediate partons are concentrated within small 
b domains. In other words, most of enhanced absorption occurs within the same parton shower, 
and is due to secondaries produced during the evolution with the same impact parameter b. On 
the other hand, a detailed analysis [14] of HERA data shows that the value of the saturation 
scale Qs{b) decreases rapidly with increasing b. {Qg is the inverse size of the dipole for which 
absorptive corrections become important.) This indicates that parton-parton correlations are 
too strong in the present model. Therefore we consider our results are close to the maximum 
possible gap suppression. 

6 Possible role of NLO effects 
6.1 Threshold factors 

The evolution equation for il^, (2), and the analogous one for Q", are written in the leading 
ln(l/a:;) approximation, without any rapidity threshold. The emitted parton, and correspond- 
ingly the next rescattering, is allowed to occur just after the previous step. On the other hand, 
it is known that a pure kinematical tmin effect suppresses the probability to produce two partons 
close to each other. Moreover, this tmin effect becomes especially important near the produc- 
tion vertex of the heavy object. It is therefore reasonable to introduce some threshold rapidity 
gap. Ay and to compute S'^^jj only allowing for absorption outside this threshold interval. The 
results of Fig. 6 were calculated for Ay = 1.5. The same rapidity interval was used to separate 
low- and high- mass diffractive dissociation in the model [1]. Below, we present the results for 
the survival factors for Ay = 0, 1.5, 2.3. For Ay = 2.3 all the NLL BFKL corrections [16] may 
be reproduced by the threshold effect [2, 3, 4, 17]. Technically we can allow for a non-zero Ay 
by just shifting yj to yj = Inl/^j - Ay. 

For the CED production of a high Et photon pair at the Tevatron with the "threshold" 
Ay = 2.3 (Ay = 0) we reproduce only one half (1/7) of the value of the full gap survival 
factor S^fi — 0.05 that was used in [18] to predict the corresponding cross section. For lighter 
CED Xc production, where the 'enhanced' absorption is stronger (that is S'enh is smaller) with 
Ay = 2.3 we now obtain'^ about factor 0.35 smaller prediction than that published in [15]. For 
the CED Higgs boson production at the LHC, the model gives S^q = 0.004, 0.09 and 0.015 for 
Ay = 0, 1.5 and 2.3 respectively. 

■^We assume that the non-pcrturbative contribution is caused by the largest-size Pomeron component Pi, 
while the perturbative contribution is due to the second component P2. 
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The existing data from the Tevatron for exclusive 77 [19], dijet [20] and Xc [21] production 
are in agreement with the predictions [18, 15], indicating that actually the value of {S'^^^^ is 
larger than obtained by the 'soft' model of [1]. So, we should regard 

{SIq) = 0.015 ±0.01 (12) 

as a conservative (lower) limit for the gap survival probability in the exclusive production of a 
Higgs boson of mass 120 GeV at the LHC Recall that this effective value should be compared 
with S"^ obtained using the exponential slope B = A GeV~^. 



6.2 Black Disc Regime (BDR) 

Note that, at relatively low scales, < 3 GeV^, the NLO gluon density extracted from the 
global parton analyses [22] does not increase with 1/x. Since the absorptive cross section 
c"abs ~ Ois/Q"^ decreases with increasing faster than the gluon density grows {xg ~ {Q'^Y' with 
7 < 0.5), the role of larger scales is even smaller. Given this phenomenological observation, we 
may assume that, for ^/s > 1 TeV, the parton density which controls the 'enhanced' absorptive 
effect, does not grow with increasing energy. Hence the S^^^ factor will depend only on the 
available rapidity interval for enhanced rescattering; that is, on the M\/ s ratio, but not on the 
initial energy. 

This observation clearly contradicts the assumption that the black disc regime (BDR) will 
have been reached at the LHC energy [23]; and that the low x gluon density will be so large 
that only on the far periphery of the proton, b > 1 fm, will there be a chance to avoid an 
additional inelastic interaction, so that the gap survives. 

Recall that the BDR of [23] was obtained using leading order (LO) gluons which grow 
steeply with 1/x. This growth is simply an artefact of the absence of a LO coefficient function, 
C^^J = 0, corresponding to ^g fusion in deep inelastic scattering (DIS), and that there is no 
1/z singularity in the LO quark-quark splitting function, Pqq{z). In order to reproduce the DIS 
data, these deficiencies are compensated by an artifically large gluon density at low x values^. 
When the NLO (and NNLO) contributions are included the DIS data are described by a flat 
(or even decreasing with 1/x) gluon at the relevant moderately low scales. With such a flat 
distribution, we expect that the value of S^^^^ for central exclusive Higgs production at the 
LHC will be larger than that for the exclusive production of a pair of high Ex photons at the 
Tevatron 

S^^^'iMn > 120 GeY) > Sj^'^"^ {^r, Et > 5 GeY) > S^T'^°^{xc). (13) 

Actually, in the preceding model of soft interactions [1], the parton density grows slowly 
with energy. After the absorptive factors e"'^^/^ are included in the evolution (2), the power 

^We thank Robert Thorne for discussfons. 
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growth of the opacity Q ~ s"^ is replaced by the logarithmic asymptotic behaviour ~ In Ins 
[24]. As a result, instead of a strong BDR suppression, the value of (S^^^^) decreases by only 
20-30% when going from the Tevatron to the LHC energy, keeping the ratio M^/y^ = 0.01 
fixed. 

In our model we do not specify the nature of the soft partons. Howver, when we introduce 
BFKL we imply the partons arc gluons, and so the screening exponents were chosen assuming 
gluons in both Pomerons. In terms of QCD, at small x and low scales, the increase of the 
parton density with decreasing x found in our model may reflect the growth of the NLO sea 
quark density^ with decreasing x. This growth is less steep than that of the LO gluons, but is 
not negligible. 

Note that the central exclusive amplitudes of interest are driven by gluon-gluon fusion: for 
example, gg H, gg ^ 77 and gg — > Xc- ^i, ^ noted above, there is a screening (NLO) 
sea-quark contribution in these amplitudes, then the effective colour factor is smaller. 

So, actually, in Nature, we expect a /arger value of 5'enh- First, because of the colour-factor 
suppression when the gluons are screened by quarks, which was not included in [1]. Secondly, 
due to a much stronger 'hot-spot' parton-parton correlation in 6-space in the model of [1] than 
that indicated by the HERA data. 

7 Comparison with exclusive Tevatron data 

This conclusion is conflrmcd by the CDF exclusive data obtained at the Tevatron. The recent 
observations of central exclusive 77 [19], dijet [20] and Xc [21] production are in overall agree- 
ment with the old KMR predictions based on the global soft model of 2000 [9]. In detail, the 
data exceed the predictions [18] for 77 production, although here the present event rate is very 
low; whereas the predictions [15], based on the KMR2000 model, are a bit above the observed 
Xc production. 

How does the present soft model [1], which includes enhanced rescattering effects change 
the predictions? Clearly, the overall gap survival probabihty S'^ will be smaller than that based 
on KMR2000. However, the difference is not large for the exclusive production, pp p + A + p, 
of a heavy system A. To be precise, for Ma — 120 GeV, the 'old' soft survival prefactors are^ 

(52^ (p2)2 ^0.0015 at the LHC (KMR2000) 

{S^){Ptf ^0.0030 at the Tevatron, 
while the present soft model [1], which includes enhanced recattering, gives 

(52)(p2)2 = 0010 at the LHC (KMR2008) 

^At NLO the gluon density in this domain is approximately flat in x. 

^The results, based on KMR2000 [9], in the form {S'^){p1f were published in [15]. 
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(-5^) (P?)^ = 0.0025 at the Tevatron, 

where all these numbers have units of (GeV)^. We see that the inclusion of enhanced rescat- 
tering does reduce the exclusive pp ^ p + A + p cross section, but well within the uncertainty. 
However, for the relatively light Xc the available rapidity interval for the enhanced suppression 
is large. In this case we find that the inchision of enhanced rescattering causes a larger differ- 
ence - reducing the exclusive cross section by about |. As a result, the new 'prediction', which 
includes 6*01111, is a bit below the CDF data on central exclusive Xc production [21], rather than 
a bit above as for our old prediction [15] based on KMR2000 [9]. 

Recall that there is no theoretical or phenomenological reason to have a strong energy 
dependence of the gap survival factor S'^. So the CDF central exclusive data measured at the 
Tevatron can be used to check and to confirm the central exclusive cross sections predicted 
at the LHC. Forthcoming higher-statistics data from the Tevatron will be valuable to further 
check our predictions for exclusive processes at the LHC. 

8 Summary 

We briefly list the main findings of our analysis. 

• Enhanced rescattering of intermediate partons occuring in an exclusive amplitude is not 
negligible. It violates soft-hard factorization, modifying the input parton distributions 
used in the QCD calculation of the bare pp ^ p + A + p amplitude, in comparison with 
those measured in deep inelastic scattering. 

• However, the suppression caused by S^^^^i numerically not large for the exclusive pro- 
duction of a heavy mass system (like a Higgs boson) at the LHC. 

• We noted two results of NLO analyses of HERA data. First, 'flat' small-x gluons are 
obtained at the relevant scales in the global fits of deep inelastic and related data. Second, 
small values of the saturation scale Qsib) are obtained for impact parameters larger than 
about 0.6 fm. These observations imply that the gluon density, which gives rise to the 
enhanced absorptive correction, is not large, and that the black disc regime is not reached 
at the LHC. 

• The energy dependence of the gap survival factor is smooth. Therefore exclusive pp — > 
p + A + p data observed at the Tevatron can be used to evaluate the suppression expected 
at the LHC. 

• The present Tevatron data imply that the predictions given here (and in related papers) 
should be considered as conservative lower limits. This offers the possibility for the 
observation of Higgs bosons in exclusive processes at the LHC. 
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